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Abstract 



Many applications of magnetic resonance are limited by rapid loss of spin coherence caused by large 
transverse relaxation rates. In nuclear magnetic resonance (NMR) of large proteins, increased relax- 
ation losses lead to poor sensitivity of experiments and increased measurement time. In this paper 
we develop broadband relaxation optimized pulse sequences (BB-CROP) which approach funda- 
mental limits of coherence transfer efficiency in the presence of very general relaxation mechanisms 
that include cross-correlated relaxation. These broadband transfer schemes use new techniques of 
chemical shift refocusing (STAR echoes) that are tailored to specific trajectories of coupled spin 
evolution. We present simulations and experimental data indicating significant enhancement in the 
sensitivity of multi-dimensional NMR experiments of large molecules by use of these methods. 
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1 Introduction 



The loss of signal due to spin relaxation 1 is a major problem in many practical applications of 
magnetic resonance. An important application is NMR spectroscopy of proteins G] ■ Multidi- 
mensional coherence transfer experiments in protein NMR are characterized by large transverse 
relaxation rates. When these relaxation rates become comparable to the spin-spin couplings, the 
efficiency of coherence transfer is considerably reduced, leading to poor sensitivity and limiting the 
size of macro molecules that can be analyzed by NMR. Recent advances have made it possible to 
significantly extend the size limit of biological macro molecules amenable to study by liquid state 
NMR [4-7] . These techniques take advantage of the phenomenon of cross-correlation or interference 
between two different relaxation mechanisms [8-13] Until recently, it was not clear if further im- 
provements can be made and what is the physical limit for the coherence transfer efficiency between 
coupled spins in the presence of cross-correlated relaxation. In our recent work, using methods 
from optimal control theory, we derived fundamental limits on the efficiency of polarization transfer 
in the presence of general relaxation mechanisms [14-16]. This established that state of the art 
experiments in NMR have the potential for significant improvement. We also provided relaxation- 
optimized pulse sequences which achieve the theoretical maximum transfer efficiency for a single 
spin pair. However, in order to apply these methods to practical NMR experiments, one needs to 
simultaneously address a family of coupled spin pairs with dispersion in their Larmor frequencies. 
In the limiting cases where cross-correlation rates are either much smaller or much larger than the 
spin-spin coupling, modifying the narrow-band relaxation optimized pulses into broadband transfer 
schemes is straight-forward by use of conventional refocusing techniques. However, in experiments, 
where both coupling and cross-correlation rates are comparable, the use of conventional refocusing 
methods for making relaxation optimized sequences broadband significantly reduces the transfer 
efficiencies as these methods eliminate either the spin-spin couplings or the cross-correlation effects. 
Finding broadband transfer schemes which can achieve the efficiency of relaxation-optimized se- 
quences required the development of specific trajectory adapted refocusing (STAR) methods, where 
refocusing is performed in a moving coordinate system attached to an optimal trajectory. In this 
paper, we present these new methods and resulting broadband relaxation-optimized polarization 
transfer experiments. 
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2 Theory 



We consider an isolated heteronuclear spin system consisting of two coupled spins 1/2, denoted I 
(e.g. 1 H) and S (e.g. 15 N). We address the problem of selective population inversion of two energy 
levels (e.g. a(3 and (3(3) as shown in Fig. 1. This is a central step in high-resolution multi-dimensional 
NMR spectroscopy [T7| and corresponds to the transfer of an initial density operator I z , representing 
polarization on spin /, to the target state 2I Z S Z , representing two-spin order. 





Figure 1: The broadband transfer of polarization I z (A) to 2I Z S Z (B) corresponds to an offset- 
independent, but transition-selective population inversion of the energy levels a(3 and (3(3 



For large molecules in the so-called spin diffusion limit |17| , where longitudinal relaxation rates are 
negligible compared to transverse relaxation rates, both the initial term (I z ) and final term (2I Z S Z ) 
of the density operator are long-lived. However, the transfer between these two states requires the 
creation of coherences which in general are subject to transverse relaxation. The two principal 
transverse relaxation mechanisms are dipole-dipole (DD) relaxation and relaxation due to chemical 
shift anisotropy (CSA) of spins / and S. The quantum mechanical equation of motion (Liouville-von 
Neumann equation) for the density operator p |17| is given by 

p = it J[-i2I z S z , p] +t: k DD [2I z S z , [2I z S z ,p\] +7T kc SA [I z , [I z ,p]] +7r k CSA [S z , [Sg,p]] 

p}} +7T kf )D/CSA [2I z S z , [S z ,p}\, (1) 

where J is the heteronuclear coupling constant. The rates fcuu, kQ SA , k^ SA represent auto- 
relaxation rates due to DD relaxation, CSA relaxation of spin / and CSA relaxation of spin S, 
respectively. The rates ^^/csa ano - ^dd/csa represent cross-correlation rates of spin / and S 
caused by interference effects between DD and CSA relaxation. These relaxation rates depend on 
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various physical parameters, such as the gyromagnetic ratios of the spins, the internuclear distance, 
the CSA tensors, the strength of the magnetic field and the correlation time of the molecular tum- 
bling |17j . Let the initial density operator p(0) = A and p(t) denote the density operator at time 
t. The maximum efficiency of transfer between A and a target operator C is defined as the largest 
possible value of Trace{C^p(i)} for any time t |2l)j (by convention operators A and C are normalized). 




Figure 2: (A) Schematic representation of the magnetization vector r% = ({I x ), (Iy)i (Iz)) an d of 
the antiphase vector r 2 = ({2I X S Z ), (2I y S z ), {2I Z S Z }) in the common frame spanned by the standard 
Cartesian unit vectors i, j, and k. The vectors l\ and I2 are the projections of ri and r 2 into the 
transverse plane and 7 is the angle between l\ and l%- (B) For the optimal CROP (cross-correlated 
relaxation optimized pulse) trajectory |14j . the units vectors e\, e 2 in the direction of r% and r 2 
are orthogonal and together with e% = e\ x e 2 define a specific moving frame along the optimal 
trajectory. (C) The pulse element Ri(t) consists of a tt rotation of spin / around e\, which leaves 
ri invariant and inverts (dashed arrow) and a 7r rotation of spin S around an arbitrary axis in 
the transverse plane, which also leaves ri invariant and brings r 2 back to its initial position (solid 
arrow). Hence, i?i neither changes the ratio 1^1/ I'll nor the angle 7 which are both constants of 
motion for the optimal CROP trajectory. 



In our recent work |14j we showed that for a single spin pair IS, the maximum efficiency rj of 
transfer between the operators I z and 2I Z S Z depends only on the scalar coupling constant J and the 
net auto-correlated and cross-correlated relaxation rates of spin /, given by k a — + k T CSA and 
k c = ^dd/csa^ respectively. Here the rates k a and k c are a factor of tt smaller than in conventional 
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definitions of the rates, e.g., k a = 1/ (7rTa), where T2 is the transverse relaxation time in the absence 
of cross-correlation effects ^] El ■ The physical limit r\ of the transfer efficiency is given by ^1] 

r,= Vl + C 2 -C, (2) 

where ( 2 = (fc„ — k 2 )/(J 2 + k 2 ). The optimal transfer scheme (CROP: cross-correlated relaxation 
optimized pulse) has two constants of motion (see Figure 2 A). If h(t) and hit) denote the two- 
dimensional vectors ((I x )(t), (I y )(t)) and ((2I x S z )(t), (2I y S z )(t)), respectively, then throughout the 
optimal transfer process the ratio IZ2I/K1I of the magnitudes of the vectors I2 and l\ should be 
maintained constant at r\. Furthermore, the angle 7 between l\ and I2 is constant throughout. 
These two constants of motion depend on the transverse relaxation rates and the coupling constants 
and can be explicitly computed |14j . These constants determine the amplitude and phase of the rf 
field at each point in time and explicit expressions for the optimal pulse sequence can be derived. 
In Fig. 3 A and B, the optimal rf amplitude and phase of a CROP sequence is shown as a function 
of time for the case k c /k a = 0.75 and k a = J. 




180°- 

1 II I II I. o° I ■ 

1 t/j-1 1 t/j-1 

Figure 3: Ideal (A, B) and approximate (C, D) implementations of an on-resonance CROP sequence 
[Tl| for k a = J and k c /k a = 0.75. Panel A shows the ideal rf amplitude A(t) — — jiB[(t) / (2ir) (where 
7/ is the gyromagnetic ratio of spins /) in units of the coupling constant J and panel C shows a 
schematic representation of an approximate CROP sequence consisting of 8 hard pulses of flip angle 
a — 21.5°. Panels B and D show the phases tp(t) of the ideal CROP sequence and its hard pulse 
approximation. 
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The transfer scheme as described assumes that the resonance frequencies of a single spin pair are 
known exactly. Therefore, the above methods cannot be directly used in spectroscopic applications 
with many spin pairs and a dispersion of Larmor frequencies. In this paper, we develop methods to 
make the above principle of relaxation optimized transfer applicable for a broad frequency range, 
making these methods suitable for spectroscopy of large proteins. A straightforward method of 
converting the smooth pulse shapes (like Fig. 3A) into a broadband transfer scheme can be realized 
by the following steps. 

a) Given the optimal amplitude A(t) and phase ip(t), of the on-resonance pulse (see Fig. 3 A and 3 
B), we can approximate the smooth pulse shape as a sequence of hard pulses with small flip angles 
ctk separated by evolution periods of duration Afc (c.f. Figs. 3 C). These are DANTE-type sequences 
(delays alternating with nutations for tailored excitation) |21|. The flip angle a& at time t is just 
rt+ js ^4^ T ^ T) w ith the phase given by tp k = ip(t) (c.f. Fig. 3 D). The delays could be chosen 
in many ways. For example, they may be all equal or can be chosen so that the flip angles ak are 
equal (c.f. Figs. 3 C). 

b) Insertion of n pulses in the center of delays to refocus the transverse components of the spins 
[22] , see Fig. 4 A-C. 



Note that this method of making relaxation optimized pulses broadband is only applicable if 
one is using either just the couplings (as in INEPT or ROPE ^H] transfer) or just the cross- 
correlation effects (as in standard CPJPT or CROP Jl] transfer for J = 0) as the transfer 
mechanism. 

For example, the relaxation-optimized pulse elements (ROPE) |15|. which only use transfer 
through couplings (special case of CROP ^1] when k c = 0) can be made broadband in a straight- 
forward way as explained above. Simultaneous 7r rotations applied to spins / and S in the middle 
of the evolution periods refocus the chemical shift evolution while retaining the coupling terms (see 
Fig. 4 B). Note however, that such a pair of ir rotations will eliminate any DD-CSA cross-correlation 
effects that might be present 7 . 

On the other hand, if J is very small or k c is close to k a (in which case transfer using cross- 
correlation effects is very efficient, c.f. Eq. |2J), it is desirable to use relaxation-optimized sequences 
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Figure 4: (A) Architecture of broadband relaxation-optimized pulse sequences, consisting of N 
periods of duration Afc (gray boxes, shown in panes B-D in more detail) and hard rf pulses with 
flip angles ak and phases (fk- (B) Chemical shift refocusing scheme preserving transfer through J 
coupling but eliminating transfer though the cross-correlated relaxation rate k c . (C) Chemical shift 
refocusing scheme preserving transfer through k c but eliminating transfer though J. (D) STAR echo 
scheme preserving both transfer through k c and though J. In B, C, and D, black bars represent 
180° rotations around an axis in the x-y plane, white bars represent 180° rotations around tilted 
axes. 

which only use cross-correlation effects for transfer (special case of CROP ^3] when J = 0). Such 
a relaxation optimized transfer is characterized by a smooth rotation I z — ► I x and vice versa 
{—2I X S Z — > 2I Z S Z ). Again such a transfer can be made broadband as explained above. In this 
case the refocusing tt pulses are applied only to spin I in the center of delays (see Fig. 4 C). By 
such pulses, cross-correlation effects are retained but coupling evolution is eliminated 7 . 

Therefore the advantage of the CROP pulse sequence (which simultaneously uses both J cou- 
plings and cross-correlation effects) would be lost in using this conventional strategy to make these 
sequences broadband. The key observation for making CROP transfer broadband is that in the 
on-resonance CROP transfer scheme, the magnetization vector 

n(t) = </»>(*) i + (I v )(t) j + </*)(*) k 

always remains perpendicular (c.f. Fig. 2 A and B) to the net antiphase vector 

ra(t) = (2I x S z )(t) i + (2I y S z ){t) j + (2I z S z )(t) k, 

where i, j, and k are the standard Cartesian unit vectors (for details see Supporting Methods). Let 
ei, e-i denote unit vectors in the direction of n and ri and let 63 = e\ x ei denote the unit normal 
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pointing out of the plane spanned by e\ and e 2 . 



Let R\(t) denote a 7r rotation of spin I around e\{t) and a simultaneous it rotation of spin S 
around an arbitrary axis in the transverse plane. Observe that Ri(t) fixes the vectors r\(t) and 
r 2 (t), see Fig. 2 C. Similarly let R 2 (t) denote a ir rotation around e 2 (t) and a simultaneous it 
rotation of spin S around an arbitrary axis in the transverse plane. i? 2 (i) inverts r\{t) and r 2 (t), i.e. 
r\{t) — > —ri(t) and r 2 {t) — > —r 2 (t). We also define i?3(<) as a 7r rotation around 63(f) which also 
results in ri(i) — ► — ri(t) and r 2 (f) — > — r 2 (t). Note that these rotations are special because they 
neither change the ratio ^ 2 |/|^i| nor the angle 7 between the transverse components l\ and l 2 . 

We now show how the rotations R\ and i? 3 can be used to produce a broadband cross-correlated 
relaxation optimized pulse (BB-CROP) sequence. Given the implementation of the on resonance 
CROP pulse (Fig. 3 A and 3 B) as a sequence of pulses and delays (Fig. 3 C and 3 D), the chemical 
shift evolution during a delay A can be refocused by the sequence (c.f. Fig. 4 D) 

A A A A 

4 3 4 1 4 3 4 

The rotations R\{t) and R 3 {t) are defined using the optimal trajectory and keep changing from one 
delay to another, as the vectors r\(t) and r 2 {t) evolve. We refer to this specific trajectory adapted 
refocusing as STAR. To analyze how this refocusing works, at time instant t consider the coordinate 
system defined by ei(i), e 2 {t) and e 3 (t) (c.f. Fig. 2 B). The unit vector along z can be written as 
aei(t) + be 2 (t) + ces(t). The chemical shift evolution generator I z can be expressed as 

I z = al ei + bl e2 + cl e3 (3) 

and the evolution for time -j under the chemical shift takes the form cxp{— i u>(al ei +bl e2 +c/ e3 )A}. 
Assuming that the R3 rotation is fast, so that there is negligible chemical shift evolution (and 
negligible relaxation) during the i?3, the sequence ^ R3 — produces the net evolution 

exp{-i uj(al ei + bl e2 + cl e , 6 )^} R 3 exp{-i u>(al ei + bl e2 + cl e3 )^} 

= R 3 exp{-i Lu(-al ei - bl e2 + cl e3 )^} cxp{-i w(al ei + bl &2 + cl e3 )^}. 

For delays A <c the effective evolution can be approximated by R3 exp{— i uj c I e3 A}. Now the 
rotation i?i can be used to refocus the remaining chemical shift evolution due to I e3 by the complete 
STAR echo sequence ^R 3 ^R\^Rz^. The effective evolution during the period A 

Ri cxp{i wc/ e3 y}exp{-i wc7 e3 y} w Ri, 



i.e. chemical shift evolution is eliminated. Note, we assume that the frame e\, e% does not evolve 
much during the four A periods so that the two R3 rotations are approximately the same. Under 
this STAR sequence, the general coupling evolution exp{— i2?:JI z S z } and the general Liouvillian 
evolution (containing cross correlation effects) is not completely preserved. Inspite of this, the 
evolution of ri(t) and ^(i) for the CROP trajectory is unaltered. This is because, for this specific 
trajectory, the magnitude of the transverse components h(t) and hit) and the angle 7 between 
them is not changed by application of these tailored refocusing pulses. Since all evolution is confined 
to transverse operators, the efficiency of the BB-CROP pulse is unaltered by application of STAR 
refocusing pulses. 



3 Practical considerations 



180° rotations around tilted axes as required by the operations and R3) of the STAR echo 
method can be realized in practice by off-resonance pulses. For example, a 180° rotation around an 
axis forming an angle 9 with the x axis can be implemented by a pulse with an rf amplitude v\ and 
offset ^ ff = v\ tan0 with a pulse duration t p — l/(2h , c g), where f^g = v\ + u^ s . At the start of the 
pulse, we assume that both the on-resonance and off-resonance rotating frames are aligned. In the 
off-resonance rotating frame the axis of rotation does not move. After the pulse, the off-resonant 
rotating frame has acquired an angle of o ff = 2ttv st relative to the on-resonance frame. For a 
pulse sequence specified in the on-resonance rotating frame, this can be taken into account by adding 
the phase </> ff acquired during a given off-resonant 180° pulse to the nominal phases of all following 
pulses on the same rf channel (In the sequence provided in supporting methods, this correction 
has been incorporated). Alternative implementations of rotations around tilted axes by composite 
on-resonance pulses would be longer and could result in larger relaxation losses during the pulses. 

Under the assumption of ideal impulsive 180° rotations (with negligible pulse duration and negligible 
rf inhomogeneity) , the STAR approach realizes a broadband transfer of polarization that achieves the 
optimal efficiency as given in |14| . However, spectrometers are limited in terms of their maximum rf 
amplitude and homogeneity of the rf field. Therefore in practice, pulses have finite widths and hence 
evolution (especially relaxation) becomes important during the pulse duration. The effect becomes 
pronounced as the number of 180° pulses is increased in order to keep the refocusing periods 
short for a better approximation to the on-resonance CROP pulse. We observe that after a point 
the loss caused due to relaxation during pulse periods overshadows the gain in efficiency one would 
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expect by finer and finer approximations of the ideal CROP trajectory. Furthermore, dcphasing due 
to rf inhomogeneity increases as the number of 180° pulses is increased. Therefore one is forced to 
find a compromise between loss due to a large number of 180° pulses versus: (a) loss of efficiency 
due to a coarser discretization of the CROP pulse, (b) reduced bandwidth of frequencies that can be 
refocused by an increased duration of the refocusing periods. When the number of refocusing periods 
becomes small, it is important to find a good way to discretize the CROP pulse so as to maximize the 
efficiency of coherence transfer that can be achieved by a pulse sequence with a prescribed number of 
evolution periods. We have developed rigorous control theoretic methods based on the principle of 
dynamic programming |26| to efficiently achieve this discretization (see Supporting Methods). This 
helps us to compute optimal approximations of CROP pulse sequences as a series of a small number 
of pulses and delays very efficiently. 



A B 




Figure 5: The buildup of antiphase vectors ri is shown for 11 different offset frequencies in the range 
of ±5 J during (A) a selective CROP (without STAR echoes) and (B) a corresponding BB-CROP 
(with STAR echoes) sequence consisting of 12 periods A&. The sequence was optimized for k a = J 
and k c /k a = 0.75 and a maximum rf amplitude of 67 J. 



Fig. 5 A shows the buildup of antiphase vectors T2 for 11 different offset frequencies in the 
range of ±5J (corresponding to ±1 kHz for J rs 200 Hz) during a CROP sequence consisting of 12 
periods A& without STAR echoes. As expected, the optimal transfer efficiency is only achieved for 
spins close to resonance. In contrast, a corresponding BB-CROP experiment with STAR refocusing 
produces efficient polarization transfer for a large range of offsets (c.f. Fig. 5B). Figure 6 shows how 
the BB-CROP sequence "locks" the angle 7 between I2 and l\ (c.f. Fig. 2) near its optimal value as 
given by on-resonance CROP pulse. 

We have carried out extensive simulations to study the loss in efficiency due to a large number 
of 180° pulses for realistic as well as hypothetical values of rf amplitudes. Fig. 7 illustrates how 
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the offset dependence of the transfer efficiency rj is effected by increasing the number of STAR echo 
periods both in the absence and presence of rf inhomogeneity. From the figures it is clear that one 
has to find an optimal number of evolution periods that gives the best performance for given system 
parameters like maximum rf amplitude, rf inhomogeneity and the bandwidth one desires to cover. 



A B 




0.4 0.8 1.2 0.4 0,8 1.2 f/j- 1 



Figure 6: Evolution of the angle 7 during the selective CROP (A) and BB-CROP (B) sequence as 
in Fig. 5 for the on-resonance case (red curves) and for an offset of -3 J. The optimal value of 7 to 
be maintained during the CROP trajectory is indicated by dashed lines. 



It is important to note that with high-resolution spectrometers, equipped with more rf power, 
relaxation losses during pulse periods can be made very small. This is illustrated in Figs. 7 A and B, 
assuming a maximum rf amplitude on the / channel of 500 J and 67 J, respectively, corresponding 
to 180° pulse durations of 5 /is and 39 fis (typical value for 13 C pulses) for the J 194 Hz coupling 
constant of the 13 C- 1 H spin pair of a model system |14M15| . (vide infra). For short 180° pulses (large 
rf amplitude) during which relaxation losses become small, a larger number of refocusing pulses has 
the largest bandwidth and approaches the ideal CROP efficiency most closely (c.f. green curve in 
Fig. 7 A). 

The refocusing sequence ~i?3— -Ri^i23~ as described in the theory section is not the only 
STAR refocusing scheme for making CROP sequences broadband. For example, -ji?3-ji?2-j-^3T 
or -ji?2-j-Rix^2^ will also perform STAR refocusing. However as indicated above, in practice it 
may be necessary to have A as large as possible, in which case one should try to refocus the largest 
of the components a, b, c of the chemical shift generator I z (c.f. Eq. more often during the 
refocusing cycle A. For example, the choice of the refocusing cycle presented in the paper is optimal 
for the values of k c /k a = 0.75 and k a / J — 1, in which case the vector 63 is mostly in the x-y plane 
and hence the magnitude of component c is smaller than the magnitude of a or b. Therefore it is 
of advantage to refocus a and b more often by performing i?3 rotations (the R3 rotation refocuses 
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the a and b components) and hence the choice of the sequence. Since there are no rotations of spin 
S during the application of R3 pulses, the total number of pulses and the resulting effects due to 
rf inhomogeneity are minimized. Dephasing losses due to rf inhomogeneity of the 180° pulses (e.g. 
applied to spin S) can be further reduced by choosing appropriate phase cycling schemes |23l 1241 125] . 




Figure 7: Offset dependence of the transfer efficiency rj for system parameters corresponding to the 
13 C- 1 H moiety of 13 C sodium formate in glycerol with k a = J and k c /k a = 0.75 ^1] (A, B) and 
corresponding to the 1 H- 15 N moiety of a protein with a a rotational correlation time of 70 ns with 
k a = 0.8 J and k c /k a = 0.73 [7] (A and B). A maximum rf amplitude on the I channel of 500 J (A), 
67 J (B), and 550 J (C) is assumed, corresponding to a hypothetical 180° ( 13 C) pulse duration of 5 
/is (A), realistic on-resonance 180° ( 13 C) pulse duration of 39 /is (B) and a 180°( 1 H) pulse duration 
of 10 us (C) for a 1 H- 15 N coupling of J = 93 Hz. Blue, red and green curves represent BB-CROP 
sequences with 4, 8 and 12 STAR echo periods respectively (for details, see Supporting Methods). 
Solid and dashed curves correspond to simulations in the absence and presence of rf inhomogeneity, 
respectively, assuming a Gaussian rf distribution with a full width at half height of 10%. 



In many cases it might also be possible to cut down relaxation losses by suitable implementation of 
the 180° pulses. For example, in the presence of a large contribution of the dipole-dipole mechanism 
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to the transverse relaxation rates, synchronization of I and S rotations can be used to create trans- 
verse bilinear operators such as I X S X which commute with I Z S Z . This way some of the losses might 
be prevented when the antiphase magnetization is passed through the transverse plane during its 
inversion by pulses. 



4 Experimental results 



In order to test the BB-CROP pulse sequence, we chose an established model system I15j . 
consisting of a small molecule ( 13 C-labeled sodium formate) dissolved in a highly viscous solvent 
(( 2 Hg) glycerol) in order to simulate the rotational correlation time of a large protein. Both the 
simplicity and sensitivity of the model system makes it possible to quantitatively compare the 
transfer efficiency of pulse sequences and to acquire detailed offset profiles in a reasonable time. 
Because of its large chemical shift anisotropy and the resulting CSA-DD cross-correlation effects, we 
use the 13 C spin of 13 C-sodium formate to represent spin / and the attached 1 H spin to represent 
spin S with a heteronuclear scalar coupling constant of J = 193.6 Hz. At a temperature of 270.6 K 
and a magnetic field of 17.6 T, the experimentally determined auto and cross-correlated relaxation 
rates of spin / were k a ~ J and k c « 0.75 k a (solvent: 100% ( 2 Hg) glycerol). For a given pulse 
sequence element, the achieved transfer efficiency of 13 C polarization I z to 2I Z S Z was measured by 
applying a hard 90° proton pulse and recording the resulting proton anti-phase signal (initial 1 H 
magnetization was dephased by applying a 90° proton pulse followed by a pulsed magnetic field 
gradient) [To] , 



Fig. 8 shows experimental on-resonance transfer efficiencies of the conventional INEPT ^2] and 
CRIPT ^15] sequences as a function of the mixing time. The figure also shows the on-resonance 
transfer efficiency of a CROP sequence consisting of four periods (without refocusing) which 
shows a gain of 65% compared to the maximum INEPT efficiency. As expected (c.f. blue curve 
in Fig. 7 B), the broadband version of this sequence (BB-CROP) with four STAR echoes has a 
reduced transfer efficiency because of relaxation losses during the additional 180° pulses, which in 
the current experiments had relatively long durations due to the relatively small rf amplitude (13 
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Figure 8: Experimental on-resonance transfer efficiencies of the CROP (open square) and corre- 
sponding BB-CROP (filled square) sequence consisting of four periods (without and with STAR 
echoes) and a total duration of 4.2 ms. For comparison, experimental on-resonance INEPT and 
CRIPT transfer efficiencies are shown as a function of the transfer time. In the experiments, spins 
/ and S correspond to 13 C and *H in 13 C-sodium formate dissolved in ( 2 Hs) glycerol. 

kHz) of the I channel ( 13 C) (the BB-CROP pulse sequence is provided in Supporting Methods). 
Additional losses are caused by dephasing due to rf inhomogeneity, which is typically larger for the 
13 C channel (where most 180° pulses are given) compared to the X H channel. The experimentally 
determined on-resonance transfer efficiency of BB-CROP is 28 % larger than the maximum INEPT 
transfer efficiency. In Fig. 9, the experimental offset profiles of the I z — » 2I Z S Z transfer efficiency of 
BB-CROP and INEPT are compared. A reasonable match is found between the experiments and 
the simulations shown in Fig. 7 B. 




Figure 9: Experimental offset dependence of the I z — > 2I Z S Z transfer efficiency for a BB-CROP 
sequence consisting of 4 periods with STAR echoes (red) and the INEPT sequence (black). The 
resulting two-spin order 2I Z S Z was converted to antiphase coherence 2I Z S X by a hard 90° (S) pulse 
and the resulting antiphase signals are shown for 11 offsets of spin I in the range of ± 500 Hz. 
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5 Conclusion 



In this paper we introduced the principle of specific trajectory adapted refocusing (STAR), which 
was used to design broadband relaxation optimized BB-CROP pulse sequence. We would like to em- 
phasize again that with increasing rf amplitudes, the efficiency of the on-resonance cross-correlated 
relaxation optimized pulse can be closely approached by the BB-CROP sequences. As future spec- 
trometers are equipped with more rf power, we can significantly reduce the duration of 180° re- 
focusing pulses, which are the major bottleneck in BB-CROP achieving the maximum efficiency. 
Based on our simulations, we expect immediate gains in NMR spectroscopy of large proteins by use 
of the proposed BB-CROP pulses. For example, in the HSQC experiment involving 1 H and 15 N, 
with maximum rf amplitudes corresponding to 12 (j,s 1 H 180° pulses and 40 fis 15 N 180° pulses, 
we expect up to 70% enhancement in sensitivity over a reasonable bandwidth compared to state 
of the art methods. With currently available rf amplitudes, in many applications it might even be 
advantageous to use broadband versions of ROPE or optimal CRIPT (special case of CROP where 
J = 0). In these cases, we only use 180° pulses in the center of each evolution period and hence loose 
less due to relaxation during the pulses (of course, as pointed out earlier, in these cases in order to 
do a broadband transfer, we will necessarily eliminate either J couplings or cross-correlation ). In 
Figs. 8 and 9 we have not compared the sensitivity of BB-CROP with CRINEPT jjj as the latter is 
not broadband for the transfer I z — » 2I Z S Z . Similar to the on-resonance CROP pulse, we have found 
that the BB-CROP pulse sequence is robust to variations in relaxation rates. Finally, the ability of 
the BB-CROP sequence to achieve the maximum possible transfer efficiency over a broad frequency 
range by use of high rf power provides a strong motivation to build high-resolution spectrometers 
with short 180° pulses. 
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6 Supporting Methods 



6.1 Orthogonality of inphase and antiphase magnetization along CROP 
trajectories 

Let ri(t),T2(t) represent the inphase and antiphase magnetization vectors as defined in the paper. 
Let l\ and I2 represent their transverse components as shown in Fig. 1 of the main text. For the 
sake of simplicity of notation, we will also use 1% and I2 to denote the magnitudes of these transverse 
vectors, as the true meaning will be clear from the context. In ^3], it was shown that the CROP 
transfer I z — » 2I Z S Z has the following properties. Throughout the transfer, the angle 7 between 
vectors l\ and I2 is constant and the ratio ^- is constant at the value 77, where r\ is the optimal 
efficiency of the CROP pulse. The two constants of motion completely determine the amplitude A 
and the phase cf>, the CROP pulse makes with the vector l\. Furthermore 7 and 77 satisfy ^1] 



Using zi(t) — (I z )(t) and Z2 (t) = (2I z S z )(t), the inner product J between r\{t) and T2(t) can be 
expressed as 



We now compute ^ along the CROP trajectory using the following identities, where the time 
dependence of the quantities, l\, I2, Zi, Z2, A, <f> is implicit. 



- cos(# — 7) + r\ cos(6* + 7) = — , 



(4) 




J = zx{t)z 2 {t) + h{t)l 2 {t) cos( 7 ). 



dzi 
~~dt 
dz 2 
~~dt 
dh 
~dt 
dh 
dt 



2ttA z\ sin cf> — ttJ(£, h ~ X h cos(# + 7)) 



2ttA l 2 sin(7 — </)) 



2ttAIi sincj) 



27rAsin(7 — (f)) z 2 — 7rJ(£ I2 — X h cos(9 — 7)) 
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Along the CROP trajectory, §^ = implies [H| 

A sin 7 (— cos</> — — cos(7 — 4>)) = 0. 
n h 

Also along the CROP trajectory |^ = rj. Using (Eq. 0J, we then obtain that along CROP trajectory 

^- cos(# - 7) + % cos(9 + 7) = — ■ 
h h X 

This then implies dJ ^ = 0. Since J(0) = as ^2(0) = 0, we obtain that throughout the transfer 
J{t) = 0. 

6.2 Dynamic Programming method for finding optimal sequence of flips 
and delays 

We now explain the method of dynamic programming |2(i| for finding the optimal sequence of 
pulses and delays that best approximates relaxation optimized pulse sequences. The method is 
best illustrated by considering the simpler case when there is no cross-correlation in the system. 
In absence of cross-correlation, the relaxation optimized transfer of I z — > 2I Z S Z is characterized by 
gradual rotation of the operator I z — » I x , followed by the rotation 2I y S z — > 2I Z S Z Let r\ 

be the magnitude of in-phase terms, i.e., r\ = (I x ) 2 + (Iz) 2 - Let /?i be the angle r\ makes with 
the transverse plane, i.e. /3i = cos -1 (see Fig. 1). Let T2 measure the magnitude of the total 
antiphase terms, i.e., r| = (2I y S z ) 2 + (2I Z S Z ) 2 and let f3 2 = cos" 1 {2I ^ z) (see Fig. 1). Using rf 
fields, we can exactly control the angle [3\ and P2 and these are thought of as control parameters 
(see Fig. 1). During the evolution of relaxation optimized pulse sequence JH] one of the (3\ or P2 is 
zero, so we assume (/?i,/3 2 ) £ ([0, f],0) U (0, [0, §]). 




Ui </ x > U 2 <2/ y 5 z > 

Figure 10: Representation of the system variables 7*1, the angles j3i, P2, and of the control 
parameters u\ — cos/3i, u 2 — cos/32 in terms of the expectation values (I x ), (Iz), {2I y S z ), and 
(2I Z S Z ). 
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Now suppose, we only have one evolution period, consisting of a pulse and delay, at our disposal. 
We can compute this optimal pulse and delay so that at the end of the evolution period, r 2 is 
maximized. Starting with r\,r 2 as defined and a given choice of @\, /3 2 and r, the values of ri(r) 
and r 2 (r) at the end of the period r are 

r\(r) = exp(— 27rfc a r)(ri cos (3\ cos(tt Jt) — r 2 cos /3 2 sm(ir Jt)) + (7*1 sin/3i) 2 (5) 
r 2 (r) = exp(— 27rfc a T)(r2 cos/32 cos(7rJr) + ri cos/3i sin(7r Jt)) 2 + (r 2 sin/3 2 ) 2 . (6) 

We write this as (t"i(t), ^(t)) = f(ri, r 2 ; fii, 02i t)- We can maximize the expression in equation 
© and find the optimal value of t,(3%,02 an d also the largest achievable value r 2 (r). This value 
depends only on the initial value r\ and r 2 and we call it V\(r\, r 2 ), the optimal return function at 
stage 1 starting from r± , r 2 . This optimal return function represents the best we can do starting 
from a given value of n and r 2 given only one evolution period. 

Given two evolution periods, then by definition V%(ri,r2) — maxp lt 2 , T V\ {f(ri, r 2 ; 01, 02, t)). 
The basic idea is, since we have computed V\ for various values of r\ and r 2 , we can use it to 
compute V2. In general then 

V n {r\,r 2 ) = max K_i(/(ri,r 2 ;^i,/3 2 ,r)). (7) 

Thus the dynamic programming proceeds backwards. We first compute the optimal return functions 
V\ followed by V 2 and so on. Computing V/c(ri,r 2 ), also involves computing the best value of the 
control parameters (3i,(3 2 ,t to choose for a given value of state ri,r 2 at stage k. We denote this 
optimal choice as 0* (r 1; r 2 , k), (3 2 ( r ii r 2ik), t* (n, r 2 ,k), indicating that the optimal control depends 
on the state ri,r 2 and the stage k. 

In practice, the algorithm is implemented by sampling the [0,1] x [0,1] square in the r 1 ,r , 2 plane 
uniformly into say 100 points. Each of these points correspond to a different value of ri,r 2 . By 
maximizing the expression in equation we can compute the optimal /3i,/3 2 ,r for each of these 
points. This would give us 0*{r\, r 2 , 1), 2 ( r i> r 2, 1), T*(r\, r 2 , 1) and also Vi(ri,r 2 ). Now to find 
V^ri, r 2 ) at these points, we sample the control space /3 2 ) £ ([0, § ], 0) U (0, [0, |]) and r € [0, jj] 
uniformly and compute the value f(r%, r 2 , fix, (3 2 , t) for all of these samples 0\, /3 2 , r and choose the 
one that has the largest value V\(f(r\, r 2 ; /3i, /3 2 , r)) .This then is the best choice of control parameters 
if there are two evolution periods to go. We also then obtain V 2 (ri,r 2 ) = Vi(f(r 1} r 2 , /?*, /? 2 , r*)). 



19 



We can continue this way and compute V r „(ri,r 2 ). Now to construct the optimal pulse sequence 
consisting of N evolution periods, we just look at the value /3*(1, 0, N), /3|(1, 0, N), r*(l, 0, N) and 
evolve the system according to these parameters and get r\ and r 2 at beginning of stage N — 1. 
But we also know /3* (ri , r 2 , ./V — 1 ) , /3 2 , r% , N — 1 ) , t* (ri , r 2 , JV — 1 ) which is then used to evolve 
the system for one more step and so on. From the sequence /3*(ri, r 2 , fc), /3 2 r 2, T *( r i, r 2, 
fc = l...iV, the optimal flip angles can be immediately determined. 



20 



6.3 BB-CROP Pulse sequence parameters 



Table 1: Parameters of a BB-CROP sequence (c.f. Fig. 4 A and D) consisting of four STAR 
echoes optimized for k a /J — 1, k c /k a — 0.75 and amplitude A — 67 J f» 13 kHz of the / channel 
for J = 193.6 Hz. The DANTE-type on-resonance pulses applied to spin / are denoted a/.. The 
Table only specifies explicitly the parameters for the spin / rf channel, however note that the pulse 
elements R\(t) require a simultaneous hard 180° rotation of spin S around an axis in the transverse 
plane. In our experiments, the phases of the four 180° (S) pulses were chosen according to the XY-4 
cycle 0°, 90°, 0°, 90° [22] in order to reduce the effects of rf inhomogeneity of the S pulses. 



tvDe 


duration [/is] 


offset [kHz] 


Phase [deg] 


a 


5.9 







Ai/4 


300.2 






Rait) 


36.1 


-4.77 


119.5 


Ai/4 


300.2 






Ri(t) 


12.7 


37.13 


20.9 


Ai/4 


300.2 






R 3 (t) 


37.6 


-2.76 


224.6 


Ai/4 


300.2 






ai 


7.9 




393.9 


A 2 /4 


220.2 








34.1 


-6.82 


182.5 


A 2 /4 


220.2 








25.1 


15.13 


36.3 


A 2 /4 


220.2 






^ 3 (t) 


35.5 


-5.39 


223.5 


A 2 /4 


220.2 






a 2 


8.6 




340.5 


A 3 /4 


219.7 






■H3(t) 


34.5 


-6.42 


157.3 


A 3 /4 


219.7 






Ri(t) 


35.2 


5.74 


353.8 


A3/4 


219.7 








34.4 


-6.48 


136.6 


A3/4 


219.7 






«3 


7.9 




217.3 


A4/4 


301.2 






R 3 (t) 


36.7 


-4.04 


58.3 


A4/4 


301.2 








38.3 


-1.19 


269.0 


A4/4 


301.2 






R 3 (t) 


36.4 


-4.49 


344.1 


A4/4 


301.2 






Q?4 


4.2 




350.3 
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